Pion Interferometry at RHIC: Probing a thermalized Quark-Gluon-Plasma? 
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We calculate the Gaussian radius parameters of the pion- 
emitting source in high energy heavy ion collisions, assum- 
ing a first order phase transition from a thermalized Quark- 
Gluon-Plasma (QGP) to a gas of hadrons. Such a model leads 
to a very long-lived dissipative hadronic rescattering phase 
which dominates the properties of the two-pion correlation 
functions. The radii are found to depend only weakly on the 
thermalization time n, the critical temperature T c (and thus 
the latent heat), and the specific entropy of the QGP. The dis- 
sipative hadronic stage enforces large variations of the pion 
emission times around the mean. Therefore, the model cal- 
culations suggest a rapid increase of R ou t / Rside as a function 
of Kt if a thermalized QGP were formed. 

Bose-Einstein correlations in multiparticle production 
processes |l| provide valuable information on the space- 
time dynamics of fundamental interactions [Q . In partic- 
ular, lattice QCD calculations predict the occurrence of a 
phase transition at high temperature, and it is hoped that 
correlations of identical pions produced in high energy 
collisions of heavy ions lead to a better understanding of 
the properties of that phase transition (for a review on 
QGP signatures, see J3|). A first order phase transition 
leads to a prolonged hadronization time as compared to 
a cross-over or ideal hadron gas with no phase transition, 
and has been related to unusually large Hanbury-Brown- 
Twiss (HBT) radii J|-^] . The phase of coexisting hadrons 
and QGP reduces the "explosivity" of the high-density 
matter before hadronization, extending the emission du- 
ration of pions This phenomenon should then de- 
pend on the hadronization (critical) temperature T c and 
the latent heat of the transition. For recent reviews on 
this topic we refer to Jt],|| . 

Here, we investigate if and how HBT radii, character- 
izing the pion source in the final state when all strong 
interactions have ceased, depend on the properties of the 
QGP and the hadronization temperature. The QGP is 
modeled as an ideal fluid undergoing hydrodynamic ex- 
pansion with a bag model equation of state, eventually 
hadronizing via a first order phase transition |^,^|. For 
simplicity, cylindrically symmetric transverse expansion 
and longitudinally boost-invariant scaling flow are as- 
sumed This approximation should be reason- 
able for central collisions at high energy, and around 
midrapidity. The model reproduces the measured px~ 
spectra and rapidity densities of a variety of hadrons at 



■/s = 17 AA GeV (CERN-SPS energy), when assuming 
the standard thermalization (proper) time Tj = 1 fm/c, 
and an entropy per net baryon ratio of sJpB = 45 [fUM^l. 
Due to the higher density at midrapidity, thermalization 
may be faster at BNL-RHIC energies - here we assume 
Ti = 0.6 fm/c and s/ps = 200. (With these initial con- 
ditions preliminary results on the multiplicity, the trans- 
verse energy, the p^-distribution of charged hadrons, and 
the vlv ratio at y/s = 130A GeV are described quite 
well po|-|l2|; below it will be shown that HBT correla- 
tions of pions at small relative momenta do not depend 
sensitively on these initial conditions.) The energy den- 
sity and baryon density are initially distributed in the 
transverse plane according to a so-called "wounded nu- 
cleon" distribution with transverse radius Rt = 6 fm. 
For further details, we refer to refs. ]lC)|-|r2"|] . 

We shall first discuss the radii of the ir~ir~ correla- 
tion functions at hadronization. From the hydrodynam- 
ical solution in the forward light cone we determine the 
hadronization hypersurface cr M . On that surface, the two- 
particle correlation function is given by 



C 2 (pi,p 2 ) = 1 
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da ■ Ke irrq f (u ■ K/T) 
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This assumes a chaotic (incoherent) and large source. / 
denotes a Bose distribution function. The normalization 
factor J\[ is given by the product of the invariant single- 
inclusive distribution of ir~ evaluated at momenta pi and 
P2, respectively, i// 1 denotes the four- velocity of the fluid 
on the hadronization surface cr M , and = (p 1 ^ +p%)/2, 
of = p\ — p% are the average four-momentum and the 
relative four-momentum of the pion pair, respectively. 
For midrapidity pions K\\ = q\\ = 0. Note that eq. (|l|) 
accounts for the direct pions only but not for decays of 
hadronic resonances (like p° — > 7r + 7r~ etc.), which are 
known to affect the correlation functions gQUJl^] and 
which will be included below. 

One usually employs a coordinate system in which the 
long axis (z) is chosen parallel to the beam axis, where 
the out direction is defined to be parallel to the transverse 
momentum vector Kj = (pit + P2t)/2 of the pair, and 
the side direction is perpendicular to both. 

From eq.(hl), the inverse widths of the correla- 



tion function are obtained as R n 



and 



R sidc = Vln2/g* dc where g* ut , g* idc are defined by 
C2(g* ut ,9 si dc = 0) = C 2 {q* sidc ,q out = 0) = 3/2. Due to 
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the definition of the out and side direction, i? ut probes 
the spatial and temporal extension of the source while 
i? s idc only probes the spatial extension. Thus the ra- 
tio i?out/-Rside gives a measure of the emission dura- 
tion (see also eqs.(3)-(5) and discussion below). It has 
been suggested that the ratio Rout/ Rsi&e should increase 
strongly once the initial entropy density Si becomes sub- 
stantially larger than that of the hadronic gas at T c Q . 
Indeed, Fig. [l] shows that R on t/Rside is much larger if 
T c is low, such that entropy conservation dictates a long 
hadronization time. The closer T c is to the initial tem- 
perature T<(« 300 MeV for the BNL-RHIC initial con- 
ditions), the faster T c is reached from above. For 1+1 
dimensional isentropic scaling expansion the time to com- 
plete hadronization is given by th/ti = Si/ sh(T c ), where 
sh(T c ) is the entropy density of the hadronic phase at T c . 
Of course, sh(T c ) increases with T c and so the hadroniza- 
tion time decreases. As -Rout grows with the duration of 
pion emission pj, it must therefore decrease as T c in- 
creases. Also, i?out/-R s ide decreases towards large Kt as 
transverse flow reduces the scale of spatial homogeneity; 
that effect is stronger for lower T c as transverse collective 
flow has more time to develop during the lifetime of the 
QGP. 



tic contribution as well as resonance excitation, which 
dominates the cross section, via 
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FIG. 1. Ratio of the inverse widths of the correlation func- 
tion in out and side direction at hadronization, as a function 
of K T ; for T c ~ 160 MeV and ~ 200 MeV, respectively. 

We now proceed to include hadronic rescattering fol- 
lowing hadronization. To describe the evolution of 
the hadrons towards freeze-out obviously requires to 
go beyond perfect-fluid dynamics. Here, we employ a 
semi-classical transport model that treats each particu- 
lar hadronic reaction channel (formation and decay of 
hadronic resonance states and 2 — > n scattering) ex- 
plicitly ]T^ , p^[ . The transition at hadronization is per- 
formed by matching the energy-momentum tensors and 
conserved currents of the hydrodynamic solution and of 
the microscopic transport model, respectively (for de- 
tails, see Q). The microscopic model propagates each 
individual hadron along a classical trajectory, and per- 
forms 2 — ► n and 1 — > m processes stochastically. E.g., 
the total meson-meson cross section includes a 5 mb elas- 
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with the total and partial -y/s-dependent decay widths 
r to t and Tji^mm' ■ (Here, -y/s refers to the cm-energy 
of the hadron- hadron scattering.) The full decay width 
r^t (to) of a resonance is defined as the sum of all par- 
tial decay widths and depends on the mass of the excited 
resonance. The sum extends over all resonance states 
which have a decay channel into MM' . The pole masses 
and partial decay widths are taken from [ jl6| . The tttt 
cross section for example is dominated by the p(770), 
with additional contributions from higher energy states. 
irK scattering is either elastic or proceeds through forma- 
tion of mainly a K * (892) resonance. Meson-baryon elas- 
tic and resonant scattering is also taken into account. In 
this way, a good description of elastic and total pion cross 
sections in vacuum is obtained |15|. Collective medium- 
induced effects on the pion scattering |l7| are neglected 
at present because in a thermalized state at T = T c there 
is less than one ir~ per phase space cell d 3 a;d 3 p/(27r) 3 . 

The distribution of freeze-out points of pions in the for- 
ward light-cone is rather broad in time pq , ^9[ . Freeze- 
out occurs in a four- dimensional region within the for- 
ward light-cone rather than on a three-dimensional "hy- 
persurface" . The single-particle distribution of pions is 
not altered very much during the dissipative hadronic 
phase because hard collisions are rare [ pdp2] | . (The pres- 
sure p is small and moreover —pdV mechanical work 
is largely compensated by +TdS entropy production.) 
There are, however, numerous soft collisions fls| ], charac- 
terizing the dissipative evolution that approaches freeze 
out. That means that the hadronic system, when start- 
ing from a state of local equilibrium at hadronization, 
disintegrates slowly rather than emitting a "flash" of pi- 
ons in a rapid decay. This is fundamentally different from 
the "explosive" hadron production from the decay of a 
classical background field at the confinement transition 
via parametric resonance p(ifl . 

The solution of the microscopic transport provides the 
classical phase space distribution of the hadrons at the 
points of their last (strong) interaction. Bose-Einstein 
correlations are introduced a posteriori by identifying the 
phase space distribution at freeze-out with the Wigner 
density of the source ||8j |l3| , f2l| , S(x, K). Corrections arise 
if the pions undergo a stage of "cascading" from the 
space-time point of their production to the point of their 
last interaction p2| . However, from the above mentioned 
model for pion rescattering we find that only ~ 15 — 20% 
of the pions in the final state freeze out after an elastic 
scattering. Rather, most pions emerge from the fragmen- 
tation of a hadronic resonance, or are emitted directly 
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from the hadronization hypersurface. Thus the source 
can be considered chaotic. A more detailed discussion 



will be given elsewhere 23 1. 

We focus here on the so-called Gaussian radius pa- 
rameters, which are obtained from a saddle-point inte- 
gration over S(x,K) |p|,p4[. The correlation functions 
themselves, the value of the intercept, and the effects of 
introducing finite momentum resolutions (as in the ex- 
perimental analysis) will be discussed in |23[| . 

The HBT-radius parameters characterizing the Gaus- 



sian ansatz are 
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with ?(Kt) = x M — (x a ')(Kt) being the space- 
time coordinates relative to the momentum depen- 
dent effective source centers. The average in (^)-(|^) 
is taken over the emission function, i.e. (f)(K) = 
J d 4 xf{x)S{x, K)/ J d^xSix, K). In the osl system = 
(/3 t ,0,A), where (3 = K/E K and E K = Vm 2 + K 2 . Be- 
low, we shall cut on midrapidity pions (f3i ~ 0), thus 
the radii are obtained in the longitudinally comoving 
frame. In the absence of x-t correlations, a large du- 
ration of emission At — \f(t 2 ) increases i? ut relative to 
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FIG. 2. HBT radii 7? ou t, Rside, J?ion g , and emission dura- 
tion At at freeze-out as a function of ATr; for central colli- 
sions at RHIC (above) and SPS (below) and T c ~ 160 MeV 
(left) and T c ~ 200 MeV (right). 

Fig. H shows the HBT radii and At. Also, we have 
varied the bag parameter B from 380 MeV/fm 3 to 
720 MeV/fm 3 , corresponding to critical temperatures of 
T c ~ 160 MeV and T c ~ 200 MeV, respectively. Note 
that within the bag model this automatically corresponds 
to a large variation of the latent heat ~ AB as well, serv- 
ing the purpose of proving the (in)sensitivity of the HBT 



radii to variations of the phase transition parameters in 
this model. We also note that changing T c implies varia- 
tion of the longitudinal and transverse flow profile on the 
hadronization hypersurface (which is the initial condi- 
tion for the subsequent hadronic rescattering stage) over 
a broad range. 

It is obvious from Fig. | that the HBT radii are rather 
similar for all cases considered. They depend only weakly 
on the specific entropy, on the critical temperature and 
latent heat for the transition, on the thermalization time 
Tj, or on the initial condition for the hadronic rescatter- 
ing stage. Thus, the properties of the QGP-phase are 
not directly reflected in the HBT radii, which are essen- 
tially determined by the large space-time volume of the 
hadronic rescattering stage. Soft hadronic scattering oc- 
curs over a long time-span after hadronization []l8|| , and 
therefore the pion mean free path increases gradually to- 
wards freeze-out p5[ . Pions are emitted over a broad 
time interval and from the entire volume. Hence, the 
emission duration is large, see Fig. @. 
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FIG. 3. Rout/Rside for RHIC initial conditions, as a func- 
tion of Kt at freeze-out (symbols) and at hadronization 
(lines) . 

In Fig. H we compare -R ut/-Rsidc at hadronization to 
that at freeze out. Clearly, up to Kt ~ 200 MeV 
-Rout/-Rsidc is independent of T c , or sh{T c ), if hadronic 
rescatterings are taken into account. Moreover, at higher 
Kt the dependence on T c is even reversed: for high 
T c the i? ut/-Rsido ratio even exceeds that for low T c . 
Higher T c speeds up hadronization but on the other hand 
prolongs the dissipative hadronic phase that dominates 
the HBT radii. This is because during the non-ideal 
hadronic expansion the scale of spatial homogeneity of 
the pion density distribution increases, as the pions fly 
away from the center, but the transverse flow can hardly 
increase to counteract (similar to a Hubble expansion 
at small pressure). Therefore, after hadronic rescatter- 
ing i? ut/-Rside does not drop towards higher Kt (in the 
range Kt < Sm^). 

Experimental data ]2(|^]] for central Pb+Pb collisions 
at y/s = 17 ' AA GeV indicate smaller HBT radii than 
seen in Fig. |^, in particular at Kt < 100 MeV. All three 
radii are < 10 fm. One should keep in mind though 
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that experimental uncertainties are largest at small Kt, 
where also Coulomb, acceptance and efficiency correc- 
tions are largest and experimental momentum resolution 
effects are crucial. At Kt > 200 MeV, the calculated 
HBT radii agree better with data. The experimental mo- 
mentum resolution or the fit procedures could in effect 
reduce the extracted radii - these rather technical issues 
will be investigated in a detailed forthcoming study Q . 
On the other hand, our calculation gives a similar in- 
dependence of the three radii as seen in the SPS data: 
strong decrease of Ri on g: weaker decrease of i? s id c and 
rather flat behavior of i? out with Kt- Most striking is 
that a new preliminary NA49 data analysis |27j] shows a 
rising i? ut/-Rsidc ratio, rather similar to our results for 
SPS energy with T c ~ 160 MeV [||. 

Nevertheless, good agreement with the data might be 
not possible to achieve for i? out and i?i on g- That would 
point at too large duration of emission, At, and too 
large spatial homogeneity scale in this model scenario 
of a first order phase transition from a thermalized QGP 
to a hadron gas phase. Note that single-particle spectra 
above pt ~ m n are described quite well ]l0|-p^|. This 
emphasizes the importance of interferometry as a cru- 
cial test of the QGP model with hadronization through a 
mixed phase, as HBT probes the small momentum trans- 
fer dynamics and the spatial homogeneity scale at freeze 
out. 

At the higher RHIC energy, the model can be tested 
at even higher Kt, where experimental resolutions and 
corrections will be less of an issue. Also, higher energy 
densities will be achieved, such that the QGP model 
should in principle be more reliable. For central colli- 
sions of Au nuclei at s/s = 130 A GeV, preliminary data 
from STAR gives R out /R si dc ^ 1.1 at small K T @. Re- 
sults from RHIC at higher Kt will test whether a long- 
lived hadronic soft-rescattering stage, associated with 
the formation and hadronization of an equilibrated QGP 
state is indeed seen in heavy ion collisions at the highest 
presently attainable energies. 
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